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Abstract 

In the framework of three-point QCD sum rules, the form factors for the semileptonic decays 
of — > B s (B*)l + ui are calculated with account for the Coulomb-like Os/v-corrections in the 
heavy quarkonium. The generalized relations due to the spin symmetry of HQET/NRQCD for 
the form factors are derived at the recoil momentum close to zero. The nonleptonic decays 
are studied using the assumption on the factorization. The B c meson lifetime is estimated by 
summing up the dominating exclusive modes in the c — ► s transition combining the current 
calculations with the previous analysis of b — ► c decays in the sum rules of QCD and NRQCD. 
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1 Introduction 

For better understanding and precise measuring the weak-action properties of heavy quarks, governed 
by the QCD forces, we need as wide as possible collection of snapshots with hadrons, containing 
the heavy quarks. Then we can provide the study of heavy quarks dynamics by testing the various 
conditions, determining the forming of bound states as well as the entering of strong interactions into 
the weak processes. So, a new lab for such investigations is a doubly heavy long-lived quarkonium B c 
recently observed by the CDF Collaboration for the first time. 

This meson is similar to the charmonium and bottomonium in the spectroscopy, since it is composed 
by two nonrelativistic heavy quarks, so that the NRQCD approach || is well justified to the system. 
The modern predictions for the mass spectra of be levels were obtained in refs. 0j in the framework 
of potential models and lattice simulations. The arrangement of excitations is close to what was 
observed in the charmonium and bottomonium. However, the feature of S c -mesons is an absence 
of annihilation into light quarks, gluons and leptons due to QCD and QED, that implies the higher 
excitations decay into the low lying levels and ground state due to the emission of photons and pion 
pairs. The measured value of B c mass yet has a large uncertainty 

M Bc = 6.40 ± 0.39 ± 0.13 GeV, 
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in agreement with the theoretical expectations. 

The production mechanism for the 5 c -meson was studied in refs. 0. The most simple picture 
takes place for the production in the e + e~-annihilation, where the universal perturbative fragmen- 
tation functions can be analytically calculated for the S-, P- and D-wave levels in the framework of 
factorization for the hard production of quarks and their soft binding into the hadron, which can be 
reliably described in the potential models. In hadron collisions, the fragmentation regime takes place 
at the transverse momenta pr > m B c) and at ~ m B c the subleading terms in 1/pt or higher twists 
have to be taken into account. This can be calculated in the framework of factorization approach by a 
careful evaluation of complete set of diagrams in the given a s -order, 0(af). The non-fragment at ional 
contributions dominate at px ~ % ||. 

The measured B c lifetime 

t[B c \ = 0.461°$ ± 0.03 ps, 

agrees with the estimates obtained in the framework of both the OPE combined with the NRQCD 
evaluation of hadronic matrix elements [[5], |6|, [7| and potential quark models, where one has to sum up 
the dominating exclusive modes to calculate the total B c width [§, || 



tope,pm[-B c ] = 0.55 ± 0.15 ps. 



The accurate measurement of B c lifetime could allow one to distinguish various parameter dependen- 
cies such as the optimal heavy quark masses, which basically determine the theoretical uncertainties 
in OPE. 

At present, the calculations of B c decays in the framework of QCD sum rules were performed in 
|TT| |T2] , [T3j . The authors of JED), |TT[] got the results, where the form factors are about 3 times less 



than the values expected in the potential quark models, and the semileptonic and hadronic widths 
of B c are one order of magnitude less than those in OPE. The reason for such the disagreement was 
pointed out in and studied in ||13|| : in the QCD sum rules for the heavy quarkonia the Coulomb- 



like corrections are significant, since they correspond to summing up the ladder diagrams, where a s /v 
is not a small parameter, as the heavy quarks move nonrelativistically, t)<l. The Coulomb rescaling 
of quark-quarkonium vertex enhances the estimates of form factors in the QCD sum rules for the 
B+ — > -?/>(?7 c )Z + z/ decays, where the initial and recoil mesons are both the heavy quarkonia. In the 
framework of NRQCD at the recoil momentum close to zero one derives the spin symmetry relations 
for the form factors of semileptonic B c decays O]. In the strict limit of V\ = t>2, where denote 



the four- velocities of initial and recoil mesons, respectively, the authors of [14] found a single relation 



between the form factors Q. In [|T3| the soft limit v\ • Vi — > 1 at V\ ^ v 2 was considered, and the 
generalized spin symmetry relations were obtained for the B c — > ip{Vc) transitions: four equations, 



including that of Moreover, the gluon condensate term was calculated in both QCD and NRQCD, 
so that it enforced a convergency of the method. 

In the present paper we calculate the B c decays due to the c — > s weak transition in the framework 
of QCD sum rules, taking into account the Coulomb-like a s /v -corrections for the heavy quarkonium 
in the initial state. In the semileptonic decays the hadronic final state is saturated by the pseudoscalar 
B s and vector B* mesons, so that we need the values of their leptonic constants entering the sum rules 
and determining the normalization of form factors. For this purpose, we reanalyze the two-point sum 
rules for the B mesons to take into account the product of quark and gluon condensates in addition 
to the previous consideration of terms with the quark and mixed condensates. We demonstrate the 



4 In refs.|l5|, [l6| the relations were studied in the framework of potential models. 
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significant role of the product term for the convergency of method and reevaluate the constants fs as 
well as fs a ■ Taking into account the dependence on the threshold energy E c of hadronic continuum 
in the bs system in both the value of fs a extracted from the two-point sum rules and the form factors 
in the three-point sum rules, we observe the stability of form factors versus E c , which indicates the 
convergency of sum rules. 

The spin symmetries of leading terms in the lagrangians of HQET JT1J for the singly heavy hadrons 
(here flW) and NRQCD § for the doubly heavy mesons (here B c ) result in the relations between the 
form factors of semileptonic B c — > B^ decays. We derive two generalized relations in the soft limit 
v i • 1)2 —>■ 1: one equation in addition to what was found previously in ref.[14|. The relations are in 
a good agreement with the sum rules calculations up to the accuracy better than 10%, that shows a 
low contribution of next-to-leading l/m^-terms. 

We perform the numerical estimates of semileptonic B c widths and use the factorization approach 
[I8| to evaluate the hadronic modes. Summing up the dominating exclusive modes, we calculate the 
lifetime of B c , which agree with the experimental data and the predictions of OPE and quark models. 
We discuss the preferable prescription for the normalization point of nonleptonic weak lagrangian for 
the charmed quark and present our optimal estimate of total B c width. We stress that in the QCD 
sum rules to the given order in ot s , the uncertainty in the values of heavy quark masses is much less 
than in OPE. This fact leads to a more definite prediction on the B c lifetime. 

The paper is organized as follows. Section 2 is devoted to the general formulation of three-point 
sum rules for the B c decays with account of Coulomb-like corrections. The analysis of two-point sum 
rules for the leptonic constant of singly heavy meson with the introduction of term allowing for the 
product of quark and gluon condensates is presented in Section 3, where we also show the convergency 
of three-point sum rules with respect to a dependence on the threshold energy of continuum in the 
heavy-light system. We estimate the form factors of semileptonic B c — > B^ decays. The relations 
between the form factors of semileptonic decays as follows from the spin symmetry of HQET and 
NRQCD are derived in Section 4 in the soft limit of zero recoil. Section 5 contains the description 
how the nonleptonic decays modes are calculated and the B c lifetime is evaluated. We discuss the 
optimal estimation of lifetimes for the heavy hadrons in Section 6. In Conclusion we summarize the 
results. 



2 Three-point sum rules for the B c meson. 



In this paper the approach of three-point QCD sum rules |TJ| is used to study the form factors of 
semileptonic and nonleptonic decay rates for the c — > s transition in decays of B c meson. From the 
two-point sum rules we extract the values for the leptonic constants of mesons in the initial and final 
states. In our consideration we use the following notations: 

f A/f 2 

(0\wi 5 q2\P( P )} = (1) 
mi + ^2 

and 

(0\qn^q2\V(p,e))=ie fM M v f v , (2) 

where P and V represent the scalar and vector mesons, and m 1; m 2 are the quark masses. 

The hadronic matrix elements for the semileptonic c — > s transition in the B c decays can be written 
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down as follows: 

(B s (p 2 )\V,\B c ( Pl )) = /+(pi+p 2 ) M + /_g M , (3) 
-(B:(p 2 )\V,\B c ( Pl )) = iF v e fivaP e*»{p 1 +p 2 )<*q' 3 , (4) 

~{B:{ V2 )\A,\B c { Vi )) = F A e; + F A (e* -pi)(pi +pa) M + F A (e* ■ Vl )q„ (5) 

where g M = (pi —p 2 )^ and e M = e^{p 2 ) is the polarization vector of B* meson. and A M are the flavour 
changing vector and axial electroweak currents. The form factors f±, Fy, F A and F A are functions of 
q 2 only. It should be noted that since the leptonic current Z M = ^7^(1 + 75)^ is transversal in the limit 
of massless leptons, the probabilities of semileptonic decays are independent of /_ and F A (the t + u t 
mode is forbidden by the energy conservation). Following the standard procedure for the evaluation 
of form factors in the framework of QCD sum rules |p0 |, we consider the three-point functions 

n M bi,p 2 ,g 2 ) = i 2 J dxdye^ x -^ ■ 

(0|T{g 1 (a;)75g2(a;),^(0),6( 2 /)75c(y)}|0), (6) 
nli A ( Pl ,p 2 ,q 2 ) = i 2 j ' dxdye^*-^ ■ 

(0|r{g 1 (x)7 M g 2 (x),J ! r' A (0),6( 2 /)75c(i/)}|0), (7) 

where qi(x)^q 2 (x) and qi(x)^ u q 2 (x) denote interpolating currents for B s and B*, correspondingly. 
J^' A are the currents and of relevance to the various cases. 
The Lorentz structures in the correlators can be written down as 



n M = n+(pi +p 2 ) M + n_^, (8) 

= ^V^^al3P 2 Pu (9) 

U A U = U^g^ + Ufp^p^ + U^p^p^ + Ufp^p^ + Ufp^p^. (10) 

The form factors f±, Fy, F A and F± are determined from the amplitudes U±, Uy, H A and U A = 
|(n^ ± n^), respectively. In (8)-(10) the scalar amplitudes 11; are the functions of kinematical 
invariants, i.e. Ilj = H^p 2 ,^, q 2 ). 

The leading QCD term is a triangle quark-loop diagram in Fig. [I], for which we can write down 
the double dispersion representation at q 2 < 

W rt (Pi,Pl^ 2 ) = -TTTXl i -T 1 — ^a!'/ 2 '^ \, dsids 2 + subtractions, (11) 



where Q = —q > 0. The integration region in (|TTJ) is determined by the condition 

X^is^s^q^X^im^s^m 2 ) ' 1 ) 

where A(xi,x 2 ,x 3 ) = (x\ + x 2 — x 3 ) 2 — Ax\X 2 . The expressions for spectral densities pf ert (si, s 2 , Q 2 ) 
are given in Appendix A. 
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Figure 1: The triangle diagram, giving the leading perturbative term in the OPE expansion of three- 
point function. 

Now let us proceed with the physical part of three-point sum rules. The connection to hadrons in 
the framework of QCD sum rules is obtained by matching the resulting QCD expressions of current 
correlators with the spectral representation, derived from a double dispersion relation at q 2 < 0. 

Ri(Pi,P2,q 2 ) = / T 1 ir77~^ ~ \ ds 1 ds 2 + subtractions. (13) 

(ZTT) J (Si — p 1 )(S 2 — p 2 ) 

Assuming that the dispersion relation flI3|) is well convergent, the physical spectral functions are gen- 
erally saturated by the ground hadronic states and a continuum starting at some effective thresholds 
s\ h and s 2 h 

pf^ s ( Sl , S2 ,Q 2 ) = P r(si,s 2 ,Q 2 )+ (14) 
9{ 8l - sf ) -9(s 2 - s\ h ) ■ P r t (s 1 ,s 2 , Q 2 ), 

where the resonance term is expressed through the product of leptonic constant and form factor for 
the transition under consideration, so that 

P r(s u s 2 ,Q 2 ) = (0|67fi( 7 „)*|fl.(fl;)) (B.{B:)\F i (Cf l )\B e ) (B c \h 5 c\0)} ■ 

(27r) 2 5(s 1 - M 2 )5(s 2 - M 2 ) + higher state contributions, (15) 

where Mi j2 denote the masses of hadrons in the initial and final states. The continuum of higher 
states is modelled by the perturbative absorptive part of IT, i.e. by pi. Then, the expressions for the 
form factors Fj can be derived by equating the representations for the three-point functions IT in flTT| ) 
and (|l~3D, which means the formulation of sum rules. 

For the heavy quarkonium be, where the relative velocity of quark movement is small, an essential 
role is taken by the Coulomb-like a s /v -corrections. They are caused by the ladder diagram, shown in 
Fig. 0. It is well known that an account for this corrections in two-point sum rules numerically leads 
to a double-triple multiplication of Born value of spectral density |19|, 21| . In our case it leads to the 
finite renormalization for p, t []13|| , so that 

P° = C Pi , (16) 
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Figure 2: The ladder diagram of the Coulomb-like interaction. 
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(17) 



where v is the relative velocity of quarks in the 6c-system, 



4m&m c 



V pl-(m b -m c ) 2 ' 
To the moment, the procedure of calculations is completely described. 



;is) 



3 Numerical results on the form factors and the semileptonic 
decay widths 

We evaluate the form factors in the scheme of spectral density moments. This scheme is not strongly 
sensitive to the value of the fee-system threshold energy. In our calculations E h c c = 1.2 GeV. The 
two-point sum rules for the B c meson with account for the Coulomb-like corrections give a^(6c) =0.45, 
which corresponds to /b c =400 MeV [21]]. The quark masses are fixed by the calculations of leptonic 
constants and fy in the same order over a s . The requirement of stability in the sum rules including 
the contributions of higher excitations, results in quite an accurate determination of masses m c = 
1.40 d= 0.03 GeV and m& = 4.60 ± 0.02 GeV, which are in a good agreement with the recent estimates 
in , where the quark masses free off a renormalon ambiguity were introduced. The values of leptonic 
constants fy linearly depend on the Coulomb-exchange a s . We find a c s (cc) ~ 0.60, a c s (bb) ~ 0.37, 
which obey the remormalization group evolution with the appropriate scale prescription, depending 
on the quarkonium contents. In this way, we can extract the above values of a c s {bc) and fs c - Note, 
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that the heavy (QiC^-quarkonia constants obey the scaling relation U 




Mjjni +m 2 ) 
\m\m2 



n 



(19) 



where n denotes the radial excitation number of nS-level, and c is independent of heavy quark flavors. 

The leptonic constant for the B s meson is extracted from the two-point sum rules. The Borel 
improved sum rules for the B meson leptonic constant |23| have the following form: 



f B M B e 



-AM' 



K 2 —C{ii) J duu 2 e- UJT + (qq)(l- 



?7T.q r 2 



vr 2 r 4 



a. 



16 



7T 



(20) 



where we use (qq) = —(0.23 • GeV)' 



rrir 



0.8 GeV 2 , (^G 2 ) = 1.77- 10~ 2 GeV 4 as the central values, 



and Mb=5.2S GeV. The K-factor is due to a s -corrections. We expect it is large, but we suppose the 
appearance of the same factor in evaluating the a s -corrections to the heavy-light vertex in the triangle 
diagram. For this factor we have the following expression 



K 



z 2 e 



-i 



E C T 



z 2 e 




In z >dz, 



(21) 



where we suppose that the scale A is equal to 1.25 GeV. The dependence of K-factor on the Borel 
parameter r and the threshold energy E c is shown in Fig. || The K-factor is not sensitive to E c 
changing in the range 1.0 -v- 1.5 GeV. 




Figure 3: The K-factor dependence on the Borel parameter r and the threshold energy E c . 

We see that NLO corrections to the leptonic constant are about 40%. Using the Pade approxima- 
tion, we find that higher orders corrections can be about 30%. So, we hold the K factor in conservative 
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limits 1.4 1.7. It is quite reasonable to suppose its cancellation in evaluating the semileptonic form 
factors due to the renormalization of heavy-light vertex in the triangle diagram. 

The contribution of quark condensate term is not sensitive to the variation of (qq) in the limits 
from -(0.23 -GeV) 3 to -(0.27- GeV) 3 (this variation corresponds to the renormalization group evolution 
and the insertion of a s -corrections to this term, so that K can be putted as the overall factor). 

In the limit of semi-local duality |23], r^Owe get the relation: A(/z) = | u (fi) (the contribu- 



tion of the quark condensate term to this equation is about ~ 15% ). We introduce the renormalization 
invariant quantities 
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flM B = K 2 -(u- al f. (22) 



For AJ^ we have AJ^ = M B — m b = 0.63 GeV, and we obtain that in the semi-local duality the 
threshold energy u^XaZ = 0-84 GeV. Neglecting the quark condensate term in the leptonic constant 
we have 

3 

Since in the three-point sum rules we use the scheme of moments and search for a stable region, in 
the general Borel scheme for f B we have to consider the stability at r / with the extended region 
of resonance contribution. We expect, that the sum rules in (^) with the redefined u ren and A ren , as 
mentioned, have a stability point at r ~ 4- 

f B M B e-^ = K^Jdco . 2 e— + <ft)(l ~^f + ^(^)), (23) 
o 

where E c is already not equal to o^XaZ- Demanding a low deviation of A from A ren = 0.63 GeV, 
we find that sum rules in Eq.(p3|) can lead to the results, which are in a good agreement with the 
semi- local duality if E c — 1.1 1.3 GeV (see Fig. |J). Then the optimal value of Borel parameter 
r = r rn ~ 6.5 GeV -1 . We write down 

f B M B e~^ = K^RE\ + (qq)(l - + ^<-G 2 », (24) 

n z 16 zoo 7r 

where R denotes the average value of e _WTm . So, we find the -Ef^-dependence of f B \JM B , whereas the 
contribution of condensate is numerically suppressed, as expected from the semi-local duality. The 
results of general Borel scheme calculations of f B ^/M B ignoring the overall i^-factor are presented in 
Fig. ffl. We observe two stability regions. The stability region at r = 2 4- 4 corresponds to that of 



considered in ||23|| . The results for the leptonic constant f B obtained from this region [23[ is about 1.5 
greater than the value obtained from the stability region at r = 6 4- 7. The second region appears only 
when we introduce the term with the product of quark and gluon condensates. The similar situation 
has been observed in the NRQCD sum rules for doubly heavy baryons |26| . The product of quark and 



gluon condensates was not taken into account in [p3|j , and therefore, the intermediate stability point 
was observed only. Fixing the optimal values of f B M B in Eq.(f24|) from Fig. ^, we can invert the sum 
rules to study the dependence of A on r, as shown in Fig. ^ where the optimal values of A agree with 
the semi-local duality and the estimate A ren = M B — Note that the intermediate stability point 



r ~ 4 exhibit a low variation of A close to 0.4 GeV, which was obtained in p6| , p3|| , and usually given 
by the potential models (see, for instance, [f2?H ). 

The calculated physical quantity should be independent of parameters in the sum rule scheme. 
However, we truncate both the operator product expansion and the perturbative series for the Wilson 
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Figure 4: fg in the semi-local duality sum rules (pale curves) and in the general Borel scheme. Solid 
lines correspond to the sum rules without condensate terms, dashed lines correspond to the accounting 
for the condensates contributions. 

A , GeV 



0.8 




2 4 6 8 10 



r, GeV 

Figure 5: The dependence of A on r with the fixed values of /bV^b, which correspond to the 
stability regions at r = 6 4- 7 GeV -1 (The notations are the same as in Fig. 
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coefficients by fixed orders in the opreator dimension and a s , respectively. This fact leads to that the 
results depend on the scheme parameters, say, the Borel variable. Moreover, the physical part of sum 
rules is modelled by the contributions of resonances and a continuum term starting at a threshold, 
that introduces the dependence on the threshold value and suggests that the stability of results can 
be improved by the terms of excited resonaces in addition to the ground state. 

Let us, at first, consider the results on the leptonic constant in the limit of semi-local duality, 
which requires the stability at r — > and corresponds to the duality region containing the ground 
state, only. The sum rules show that the condensate contributions are given by the polinomials over 
r, and the leading correction at r —>■ is the term with the quark condensate. Then, we expect that 
the region of stability will extend at r > if we will add the higher condensates with the appropriate 
ratios of their values. The central values of condensates as mentioned above correspond to the results 
shown in Fig. |j. The stability of semi- local sum rules can be improved by a variation of A and ujq, 
which is not important for the current discussion. In order to clarify this statement we present the 
results of semi-local sum rules for the leptonic constant of B meson in Fig. || where we put A = 0.54 
GeV. Note, that the value of leptonic constant is the same as we have found in the pertubative limit 
of semi-local duality. Then, we state that the value of leptonic constant obtained at r — > agrees with 
the value corresponding to the stability at r « 7 GeV -1 , i.e. in the second point of local extremum 
in Fig. |4j. In order to confirm, we present also the result for the other ratio of condensate values (the 
lower value of mixed condensate and the upper value of gluon condensate in the regions mentioned 
below) in Fig. [7|. We see that the semi-local duality is completely broken at unappropriate choice of 
condensate values. 

f B VM^, GeV 3/2 



0.5 
0.4 
0.3 



0.2 
0.1 



2 4 6 8 10 

r, GeV" 1 

Figure 6: J'b in the semi-local duality sum rules (dashed curve) and in the general Borel scheme (solid 
line) with the corrected value of A, which improves the stability of result obtained in the semi-local 
duality. 

Second, we study the general Borel sum rules in the scheme, where we do not require the stability at 
t — ► 0, which means that we extend the duality region to incorporate possible excitation contributions. 
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Figure 7: fs in the semi-local duality sum rules (dashed curve) and in the general Borel scheme 
(solid line) with the condensate values, which destroy the semi-local duality. 
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Figure 8: The dependence of A on r with the fixed values of JbVMb, which correspond to the stability 
regions at r = 6 4- 7 GeV -1 in Fig. ||. 



11 



Then, we expect that the estimate of leptonic constant should result in the same value obtained in 
the semi-local duality. We see that this can be reached at the same values of condensates involved 
as well as in the same region of r, which is again given by the second extremum (see Figs. |] and ^j. 
Certainly, the extension of duality region leads to a failure of stability at low r, which is caused by 
additional terms given by the excitations beyond the control of the method. 

At higher values of r the stability would be reached by introduction of higher condensates con- 
tributing as the polinomials of higher powers. Of course, we can achieve the better stability in the 
general Borel scheme by adjusting the condensate values, but this will destroy the semi-local duality, 
that indicates the divergency of the method, while the convergency demands an appropriate choice of 
condensate values as we put. 

The criterion on the covergency of both the semi-local duality and general Borel sum rules was 
ignored in | 23| , where the greater value of leptonic constant was obtained (see Fig. |?]). 

The same notes can be done in the discussion on A. For the sake of comparison, we present the 
results in Figs. and || corresponding to the inverted sum rules for the leptonic constant given in 
Figs. [| and || 

The physical meaning of gluon condensate as it stands in the operator product expansion taken 
between the observed states is independent of the scheme of calculations. In the sum rules the 
gluon condensate contributes in the region, where the excited states are suppressed, while the ground 
state corresponding to the current under consideration dominates. So, the condensate essentially 
determines the binding energy of quarks in the hadron containing the heavy quark. The definition of 
sum rules scheme includes the parameter determining the threshold energy of continuum contribution 
in addition to the resonance term. So, the difference between the schemes of semi-local duality and 
usual Borel representation is due to the variation of duality region. So, the semi-local duality means 
the duality for the region contaning the ground state only, while the usual Borel (or moment) scheme 
explore the extended region containing several hadronic states: the ground state and its excitations. 
However, the gluon condensate essentially contributes at the scheme parameters, when the excitations 
are suppressed. Thus, its value is the same for both schemes used, i.e. for the semi-local duality 
and in the Borel scheme. The value of gluon condensate has been varied in the range (^-G 2 ) = 
(1.5 -T- 2) • 10 -2 GeV 4 . As for the mixed condesate we have used the range m 2 , = (0.72 4- 0.88) GeV 2 . 
The variation of these parameters does not change the qualitative picture for the lepton constant as 
it has been discussed, while the numerical uncertainty of its value is less than 7%. 

Numerically, multiplying the result taken from Fig. || by the K-factor we find the value fs = 
140 -T- 170 MeV, which is in a good agreement with the recent lattice results PB1 and the estimates in 



the QCD SR by other authors [p9f| . So, we can conclude that the l/m 6 -corrections are not valuable for 
fs- The uncertainty of estimates is basically connected with the higher orders in a s . For the vector 
B* meson constant f B * we put ^ = 1.11 (see |30|, |9f). 

For the leptonic constant of B s meson we explore the following relation ^ = 1.16, which expresses 
the flavor SU(3)-symmetry violation for B mesons |25 |. 

Remember, in sum rules the heavy quark masses are fixed by the two-point sum rules for bottomo- 
nia and charmonia with the precision of 20 MeV. In our consideration the quark masses are equal to 
mfe=4.6 GeV, m c =lA GeV, and we use m s =0.15 GeV, which agrees with the various estimates |31[ 



The uncertainties in the values of form factors are basically determined by the variation of 6-quark 
mass, while changing the other quark masses in the ranges m s = 0.14-=-0.16 GeV and m c = 1.35 -hi. 45 
GeV, results in the uncertainty less than 2%. In Figs. |^, and [ll] we present the results in the scheme 
of spectral density moments. 
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Figure 9: F A in the scheme of the spectral density moments; n is the number of moment with respect 
to the square of momentum in the be channel. The number of moment with respect to the square of 
momentum in the bs channel is equal to 1 at E c — 1.2 GeV. 
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Figure 10: / + in the scheme of the spectral density moments; n is the number of moment with respect 
to the square of momentum in the be channel. The number of moment with respect to the square of 
momentum in the bs channel is equal to 1 at E c — 1.2 GeV. 
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F A ,GeV m 6 




Figure 11: Fq in the scheme of the spectral density moments; n is the number of moment with 
respect to the square of momentum in the be channel, m is the number of moment with respect to the 
square of momentum in the bs channel at E c = 1.2 GeV. 

We have investigated the dependence of form factors on the bs threshold energy of continuum in 
the range E c = 1.1 4- 1.3 GeV. The characteristic forms of this dependence are shown in Figs. |T^| 
and [13|. We see that the optimal choice for the bs system threshold energy is 1.2 GeV. In Table [1] we 
present the results of sum rules for the form factors in comparison with estimates in the framework of 
potential models || [52). We see a good agreement of estimates in the QCD sum rules with the values 
in the quark model. For the sake of completeness the quark model expressions for the form factors 
are given in Appendix B. 



Method 


/+ 


f- 


F^GeV" 1 


Fq, GeV 


Ff , GeV" 1 


F*, GeV" 1 


This paper 


1.3 


-5.8 


1.1 


8.1 


0.2 


1.8 


Potential model || 


1.1 


-5.9 


1.1 


8.2 


0.3 


1.4 



Table 1: The form factors of B c decay modes into the B s and B* mesons at q 2 = 0. 

In [ffDJ the form factors were derived using the similar SR technique but without the Coulomb-like 
corrections in the be system, which enhance the form factors about three times, as we have found. 

Let us discuss the uncertainties in the sum rules and other approaches. So, the potential models 
|| [T5| with similar choices of parameters result in the form factor values, which are slightly model- 
dependent. The corresponding accuracy is about 10%. Then, we expect that the potential models 
give good reference points for the appropriate numerical values. 

The accuracy of sum rules under consideration is basically determined by the variation of heavy 
quark masses. Indeed, the significant a s correction to the leptonic constant of B s meson should cancel 
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Figure 12: The dependence of F A on the bs threshold energy E c , determining the region of resonance 
contribution. 
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Figure 13: The dependence of /+ on the 6s threshold energy E C) determining the region of resonance 
contribution. 
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the same factor for the renormalization of quark-meson vertex in the triangle diagram. The dependence 
on the choice of threshold energy in the Es-channel can be optimized and, hence, minimized. The 
variation of threshold energy in the fee-channel give the error less than 1%. The effective coulomb 
constant is fixed from the two-point sum rules for the heavy quarkonium, and its variation is less than 
2%, which gives the same uncertainty for the form factors. The heavy quark masses are determined 
by the two-point sum rules for the heavy quarkonia, too. However, their variations result in the most 
essential uncertainty. Summing up all of mentioned variations we estimate Sf/f ~ 5%. 

The enhancment of form factors in the decays of heavy quarkonium considered in the framework 
of sum rules was also found in the decays B c — > iplv [|12|, [13]]. It is important to stress that the 
current consideration removes the contradiction between the estimates in OPE, potential models and 
the values in the SR. We have calculated the form factors in the SR consistent with the values in the 
OPE and potential models. 

We use the three-point sum rules to determine the dependence on q 2 of form factors in vicinity of 
q 2 = 0, where the method works even in the approximation by a bare quark loop [pOj] . Naively, we 
expect a simple pole form 

M) = (25) 

pole 

so that the first derivative /'(0) = f(0)jp — can be evaluated in the framework of sum rules to 

pole 

estimate the "pole" mass M po i e , which may deviate from the value of physical mass of corresponding 
bound state. Of course, the bare quark loop approximation cannot be justified at q 2 ~ m 2 ole , while 
the three-point sum rules in the form of double dispersion relation cannot be explored in the region 
of resonances at q 2 > because of so-called "non-Landau" singularities indicating the presence of 
strongly bound states in the g 2 -channel. That is why we calculate the form factors as well as their 
Borel transforms at q 2 < 0. 

Numerically, we have found M po ; e =1.3-t-1.4 GeV for the form factors with the B* decay modes, 
and Mp i e — 1.8-^-1.9 GeV for the decay form factors with the B s modes. 

The semileptonic widths are presented in Table £| We have supposed the quark mixing matrix 



element \V CS \ =0.975 |]3~T]| . The mesons masses are equal to M Bc =Q.25 GeV, M Ba =5.37 GeV, M B *=hAl 
GeV |1. 



mode 


T, 10" 14 GeV 


BR, % 


B s e + u e 


5.8 


4.0 


B* s e + v e 


7.2 


5.0 



Table 2: The widths of semileptonic B c decay modes and the branching fractions calculated at r Bc = 
0.46 ps. 

These results agree with the values obtained in the framework of covariant quark model ||: 
T(B s e + v e )=4.7 ■ 10~ 14 GeV, T(B*e + is e )=7A ■ 10" 14 GeV, as we could expect looking at Table |. 
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4 The symmetry relations 



At the recoil momentum close to zero, the heavy quarks in both the initial and final states have small 
relative velocities inside the hadrons, so that the dynamics of heavy quarks is essentially nonrelativistic. 
This allows us to use the combined NRQCD/HQET approximation in the study of mesonic form 
factors. The expansion in the small relative velocities to the leading order leads to various relations 
between the different form factors. Solving these relations results in the introduction of an universal 
form factor (an analogue of the Isgur-Wise function) at q 2 — > q^ax- 
We consider the soft limit 



w 



v x -v 2 ^> 1, 



(26) 




are the four-velocities of heavy mesons in the initial and final states. The 



where < )2 = p r 12 / 

study of region (PE|) is reasonable enough, because in the rest frame of B c meson (pj* = (yjpi, 0)), the 
four- velocities differ only by a small value \p 2 \ {p 2 = (E 2 ,p 2 ), whereas their scalar product w deviates 

1 _|_ lJi>kl£ 



from unity only due to a term, proportional to the square of \p 2 \: w 



l>2 



Thus, 



in the linear approximation at \p 2 \ 

to note, that ([26]) generalizes the investigation of WM , where the case of v% 



0, relations ( p6|) are valid and take place. Here we would like 

v 2 was considered. This 



condition severely restricts the relations of spin symmetry for the form factors and, as a consequence, 
it provides a single connection between the form factors. In the soft limit of zero recoil we find 



V3 



for the four-velocity of spectator 6-quark, and 



vl + 



2mi 



Vl - v 2 



for the decaying c-quark. The matrix element of J M 
has the form 



Qi^ixq2 with the spin structure 



^(1 + ^)^(1 + ^7, 



(27) 

(28) 
{7^,757/,} 

(29) 



where Ti determines the spin state in the heavy meson Q1Q3 (in our case it is pseudoscalar, so that 
Ti — 75), r2 determines the spin wave function of quarkonium in the final state T 2 = {75,^^7^} for 
the pseudoscalar and vector states, respectively (H = P,V). The 'propagator of the light quark'f] is 
taken in a general form 



Plight = 1 + B(ji 2 - fa) + C(j) 2 + fa) + Dfafa, 



(30) 

where B,C,D are the functions of w. The quantity h in at w — * 1 is an universal factor 
independent of the spin state of meson. So, for the form factors, discussed in our paper, we have 



(P Ql Qs\Ql^Qs\P g2 Q 3 ) 

(P Ql Q 3 \Qi^Qs\v q2 Q 3 ) 



(cf • < + cl 



v$) ■ h, 



icy ■ ^ Vap e u Vi a v 2 /3 ■ h, 
(ce • e* 4 + ci • <(e • vi) + c 2 • v%{e ■ v x )) ■ h, 



(31) 
(32) 
(33) 



'More strictly, we determine the spin structure of matrix element and the term given by light degrees of freedom. 
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where 



c t = -2, 

cv = -1-B-^-, (34) 
Ami 



°1 ■ 

= 1 + 5 



2mi 
m 3 

2mi 



and i? = The rest coefficients ci^ depend on the C and D parameters. We have the symmetry 

relations for the following form factors^]: 

f+(c£-M2-c£M 1 )-f-(cZ-M2 + <£-Mi) = 0, 

F A -c v -2c £ -F v M 1 M 2 = 0, (35) 
Fq c i + c t ■ JAi{f+ + /-) = 0, 

where M.\ = mi + m 3 , A^2 = ^2 + m 3- Equating the second relation in (|35]), for example, we obtain 



2mi + 7713 4ms(mi + 7713)^ 
~2r7^ + ^ 



B = 1J'"3 + ^»3v»l -r»^ y „ 1Q 0> (36) 



where all form factors are taken at q^ax- Substituting B in first and third relations, we get / + ~ 2.0 
and /_ ~ —8.3. These values have to be compared with the corresponding form factors obtained in 
the QCD sum rules: f+iq^ax) = 1-8 an d f-iVmax) = ~ 8.1, where we suppose the pole like behaviour 
of form factors (see Eq. fl25|) ) . Thus, we find that in the QCD sum rules, relations ( p5|) are valid with 
the accuracy better than 10% at q 2 = q^ax- The deviation could increase at q 2 < g^ aa , because of 
variations in the pole masses governing the evolution of form factors. However, in B~£ — > B^l + u 
decays the phase space is restricted, so that the changes of form factors are about 50%, while their 
ratios develop more slowly. 



5 Nonleptonic decays and the lifetime 

The hadronic decay widths can be obtained on the basis of assumption on the factorization for the 
weak transition between the quarkonia and the final two-body hadronic states. For the dominant 
nonleptonic decay modes B+ — > B^it + (p + ) the effective Hamiltonian can be written down as 

Heff = ^=V cs V: d {C + (v)0 + + C»CL}, (37) 

where 

0± = (ui7„(l - 7 5 )d i )(s/y l '(l - 75)9) ± {ucf v {\ - "f 5 )dj){siY{l ~ 75)9), (38) 

where i,j run over the colors. The factors C±(ji) account for the strong corrections to the corre- 
sponding four-fermion operators caused by hard gluons. The review on the evaluation of C±(/x) can 

6 To remove an error in |l3| the analogous second relation for the B c — * J/^(j] c ) transition should have the missed 
factor 2 in front of c £ . 
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be found in ||34 | . In the present paper, dealing with the QCD sum rules in the leading order over a s , 
we explore the C± (/^-evolution to the leading log accuracy. The B+ — > B s n + amplitude, for example, 
takes the form 

A(B+ -> B s n + ) = ^V^aM^u^l - ^^{B^ (I - ls)c\B c ), (39) 

where a^p) = ^-(C + (p)(N c + 1) + C-(fi)(N c — 1)) at N c — 3 being the number of colors. In our 
calculations we put the following light meson parameters: m 7r +=0.14 GeV, m p +=0.77 GeV, ^+=0.132 
GeV, f p + =0.208 GeV. The results are collected in Table 0. 



mode 


T, 10- 14 GeV 


BR, % 


B S TT + 


15.8 aj 


17.5 


B sP + 


6.7 a\ 


7.4 


B*ir+ 


6.2 a\ 


6.9 


B* sP + 


20.0 a\ 


22.2 



Table 3: The widths of dominant nonleptonic B c decay modes due to c — >• s transition and the 
branching fractions calculated at r Bc = 0.46 ps. We put ai=1.26. 

It is worth noting that the sum of widths for transitions B+ — > B s (B*)ir + (p + ) is 10% larger than 
the width for the transition £>+ — > B S (B*) + light hadrons, which is calculated using the simple 
formula 

F[B+ -»■ B S (B* S ) + light hadrons] = N c a\(p) T[B+ -> S,(S^)e + i/ e ], 

where we neglect the contributions given by the modes with the factor a2 instead of a\. In addition, 
the deviation between these estimates can be caused by the corresponding 'bag' factor appearing in 
the formulation of factorization approach and vacuum saturation in the connection of leptonic form 
factors to the hadronic ones. The modern lattice estimates show that the 'bag' parameters are about 
7% less than 1 PSfl. 



In the parton approximation we could expect 

T[B+ -> SW + Ught hadrons] = (2C^(/i) + C 2 (fi))T[B+ -> B&e + u e ], 

which results in the estimate very close to the value obtained as the sum of exclusive modes at 
p > 0.9 GeV. The deviation between these two estimates slightly increase at ^ < p < 0.9 GeV. 
Concerning the comparison of hadronic width summing up the exclusive decay modes with the estimate 
based on the quark-hadron duality, we insist that the deviation between these two estimates is less that 
10% and, hence, it cannot be treated as an essential argument against the validity of our calculations. 

We estimate the lifetime using the fact that the dominant modes of the B c meson decays are the 
c — ► s, b — > c transitions with the By' and J/V, Vc final states respectively, and the electroweak 
annihilation []. 

We stress that in the B c — ► B s decays caused by the weak decays of charmed quark the possible 
hadronic final states are the charged mesons tt, p and K or the multi-particle states like 7r7r, tttttt 

7 The b — ► ccs transition is negligibly small in the B c decays because of destructive Pauli interference for the charmed 
quark in the initial state and the product of decay || . 
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or Kit. First, the states with the kaon are suppressed by the Cabibbo angle, and we neglect their 
contributions in the total nonleptonic width of B c (the corresponding error of estimates is about 
4%). The method for the calculation of multi-particle branching fractions was offered by Bjorken in 
his pioneering paper devoted to the decays of hadrons containing heavy quarks |35| . He supposed a 
simple relation for the yields of pions, as given by the Poisson distribution with the average value of 
pions determined by the energy release. The Bjorken's model of multi-particle yields in the decays 
does not distinguish the resonant and continuum final states as well as B s and B*. So, we check that 
the ratio of R 21T = T(B+ -> B^*h + ) /T(B+ -> B^p + ) « 0.82 calculated in the framework of sum rules 
is close to the estimate R 27r = T(B+ -> B^ ) tt + )/T(B+ -> B^tt + tt°) » 0.85 given by Bjorken. Then, 
we see that the non-resonant multi-particle states are suppressed in comparison with the resonance 
yields. The same fact can be found, once we consider the Ktt, Kwtt and K*ir, Kp branching fractions 
in the decays of D mesons as measured experimentally. This consideration confirms us that we take 
into account all of significant nonleptonic decay modes of B c . In order to estimate the contribution 
of non-resonant 3n modes of B c decays into B^ we use the Bjorken's technique, i.e. the Poisson 
distribution with the average value corrected to agree with the non-resonant 37r-modes in the decays 
of D mesons. We have chosen the following branching ratios: BR(D + — > K~tt + tt + ) = 9.0 ± 0.6% 
and BR(D + — > X~7r + 7r + 7r°| non _ resonant ) = 1.2 ± 0.6%. So, we have found n ~ ~, which means that 
BR(5+ -> SW(37r) + ) « 0.2%, while BR(B+ -> J B(* ) ( 27r ) + lnon-re S onant) ~ 3%. We see that the 
neglected modes contribute to the total width of B c as a small fraction in the limits of uncertainty 
envolved. 
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Figure 14: The dependence of B c meson lifetime on the scale \x in the effective Hamiltonian ([37]). The 
shaded region shows the uncertainty of estimates, the dark shaded region is the preferable choice as 
given by the lifetimes of charmed mesons. The dots represent the values in the OPE approach. 



The width of beauty decay in the sum rules was derived using the similar methods in fll3| : T(£>+ — >■ 
cc + X) = (28 ± 5) • 10~ 14 GeV. The width of the electroweak annihilation is taken from as 
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12 ■ 1(T 14 GeV. 

In Fig. [L4| we present the B c meson lifetime calculated in the QCD SR under consideration. We 



also show the results of the lifetime evaluation in the framework of Operator Product Expansion in 
NRQCD|,0. 

In contrast to OPE, where the basic uncertainty is given by the variation of heavy quark masses, 
these parameters are fixed by the two-point sum rules for bottomonia and charmonia, so that the 
accuracy of SR calculations for the total width of B c is determined by the choice of scale [i for the 
hadronic weak lagrangian in decays of charmed quark. We show this dependence in Fig. [14|, where 

< fi < m c and the dark shaded region corresponds to the scales preferred by data on the charmed 
meson lifetimes. The discussion on the optimal choice of scale in hadronic decays is addressed in the 
next section. 



6 Discussion on the lifetimes of heavy hadrons 

At present the ordinary prescription for the normalization point of lagrangian generating the non- 
leptonic decays of heavy quark Q, is fi ~ uiq. The motivation is the following: the characteristic 
scale is determined by the energy release given by the heavy quark mass. Therefore, we can argue 
the operator product expansion in the inverse powers of rriQ, wherein we can factorize the Wilson 
coefficients taken in the perturbative QCD and the matrix elements of operators over the hadronic 
states with \i usually posed to rriQ. This prescription is in a qualitative agreement with the current 
data on the measured lifetimes of charmed and beauty hadrons and their branching ratios for the 
semileptonic decay modes, say. 

Let us consider this issue in more details. To the moment, the analysis of decays in the QCD sum 
rules is restricted by the leading order (LO) in a s (except the Coulomb-like corrections in the heavy 
quarkonia). The corresponding parameters, the heavy quark masses, are also fixed to the same order 
(they have to be reevaluated in next-to-leading order). Therefore, for the sake of consistency, we use 
the LO expressions, which are given by the partonic approximation improved by 1/ mg-corrections in 
HQET or NRQCD. In this way, we can write down the following formulae: 

1. The semileptonic branching fraction of D° meson is given by 

BR -' |fl01 = 2 + (40) 

2. The difference between the total widths of charmed mesons is determined by the Pauli inter- 
ference in decays of D + , so that |36|] 



T[D + ]-T[D°] = cos 4 9 C G 



2 m lfD 

F Stt 



(Cj(Ai) - C 2 _^))B + 1(C» + C 2 _ {n))B 



(41) 



where 9 C is the Cabibbo angle, fo is the leptonic constant of D meson, and the 'bag' constants are 
defined by 

(D\(cT,q)(qT,c)\D) = f 2 D M 2 D B, (42) 
(D\(cT,c)(qT,q)\D) = \f D M 2 D B, (43) 

with = 7 M (1 - 7 5 ). 
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The similar expression can be derived for the beauty mesons p6 

r[B°] 



r[s H 



\v h 



be 



2 n 2 m bfB 

<-JT7 Z 



■ill 



1 



(CJ fci) - C 2 fr))B + -(Cl(fi) + C 2 _^))B 



3. The semileptonic branching fraction of B° meson is given by 

i 

BR S ; [B 



2 + 0.22 + \2C%{n) + C 2 _ (//)](! + fc) 



(44) 



(45) 



where the fraction of 0.22 is due to the ru T - contribution, while the value of k denotes the fraction 
of b — > ccs transition in the nonleptonic decays. In the same way, the average yield of charm in the 
decays of beauty mesons is equal to 



n r 



2 + 0.22 + [2C 2 + (fi) + (//)](! + 2k) 
2 + 0.22 + \2C\ (ft) + C 2 _(fx)](l + k) ' 



(46) 



As for numerical applications, one usually puts 
Hd = m c in decays of D mesons, 
[iB = m>b m decays of B mesons, 
f D m f B m 200 MeV, and 

B = B = 1 naively motivated by nonrelativistic potential models. 

At k — 0.4 ||37|| , this set (marked as SETMQ column) results in the estimates shown in Table ^ in 
comparison with the experimental data 



quantity 


exp. 


SETMQ 


SETMQ 2 


SETH 


BR sl [D% % 


8.1 ± 1.1 


15.4 


15.4 


8.6 


T[D+]-T[D°], ps- 1 


-1.56 ±0.03 


-1.26 


-0.19 


-1.53 


r[D°]/T[D a ] 


1.12 ±0.05 


1.00 


1.00 


1.11 


BR sl [B% % 


10.45 ±0.21 


14 


14 


10.2 


F[B+] -T[B% ps- 1 


-0.043 ±0.017 


-0.022 


0.024 


-0.044 


T[B°]/T[B S ] 


1.00 ±0.05 


1.00 


1.00 


1.05 


n c 


1.12 ±0.05 


1.20 


1.20 


1.12 


r[5°]/r[A 6 ] 


0.81 ±0.05 


1.00 


1.00 


0.81 



Table 4: The comparison of theoretical estimates at various sets of parameters with the experimental 
data on the decays of heavy mesons. 

First, we note that the semileptonic width of D° is well described to the given order and at chosen 
value of m c , while its branching ratio is in a valuable contradiction with the data indicating a more 
higher enforcement of nonleptonic modes. Second, the qualitative agreement of predictions with the 
measured differences of r[D + ] — r[D°] and — T[i? ] is mainly based on the assumption of B w 1. 

Recent consideration of charmed baryon lifetimes by M.Voloshin |38[ clearly drawn a conclusion that 
the naive picture of color structure as given by the potential models (i.e. the purely antisymmetric 
color-composition of valence flavors) is significantly broken. A similar statement was obtained in the 
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description of D* meson production at HERA, where the authors of |39[ found that the 0(a em a s )- 



calculations for the differential cross sections of cg-pair composing the meson, are able to reproduce 
the measured spectra, if we introduce the valuable contribution by the color-octet state in addition 
to the singlet one. So, the four-quark singlet-operator results in 

( i) = (D*\(ry»q)(q^c)\D*) • (-g, v + P -^\ Jj^, (47) 

which is reduced to 0(\) = |\I/(i)(0)| 2 in the framework of nonrelativistic potential model, where 



with e a denoting the polarization vector and \I/(i)(g) being the wave function. 
The term of color-octet was parameterized by 

(8) = (D*\(c\^q)(q\^c)\D*) ■ (-g, v + -L, (49) 



which is, in a similar way, can be represented as 0(8) = |^/(s)(0)| 2 , if we introduce the additional Fock 
state 

where the c and q fields represent the rapid valence quarks, and n a is a random color- vector determined 
by soft degrees of freedom inside the meson (i.e. by the quark-gluon sea). We have (n a rib) = S a b in 
the production, while (n a rib) = |<5 a ft in decays. 

The data on the D* production in DIS give 0(s)/0(i) — 1.3. We can analogously expect that in 
the D meson the ratio of four-quark matrix elements is close to that in D*. So, 0(u\[D\/Om[D\ ~ 1, 
where 0(i j8 )[-D] can be obtained from the above expressions for D* by the substitution of 757^ for 7^ 
and removing the transverse projector. Then we straightforwardly find 

Putting B = 2 into the SETMQ we get the values given in the SETMQ 2 column in Table f|. So, the 
choice of [i = ttlq is in a deep contradiction with the observed differences of total widths for the heavy 
mesons at the most reasonable value of B = 2. 

In this position we argue the following: There are other physical scales in the problem, which are 
characteristic for two hadronic systems in the decay process. The first system is the decaying hadron. 
The second is the transition current c — > s or b — > c, where the form factor behaviour versus the trans- 
ferred momentum is determined by the cs and be states. Those hadronic systems have the following 
scales, being the average squares of heavy quark momentum (p 2 ), which are phenomenologically equal 
to 

/4 = /4 = 2TA, (52) 

where according to the potential models A ~ 0.4 GeV is the binding energy of heavy quark (i.e. 
the constituent mass of light quark), T ~ 0.45 GeV is the average kinetic energy in the system. T 
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determines the 2S-1S splitting and it is approximately independent of quark flavors. Analogously, we 
put 

2TA S , (53) 



2 _ 2 
ftcs ~ ^bs 



where A s = A + (tod s — mo) = A + {ms a — ^lb) — 0.5 GeV. For the b — > c current we put 

H% = 2Tm bc , 

with rribc — m&ra c / {m h + m c ) ~ m B m£i/ {m B + m D ) « 1.3 GeV. 

Let us suppose that the decay scale is given by the following combinations: 



(54) 



As 



2 

f^bs ■ Hcbi 



(55) 



At B = 2, this set of parameters with fB — fo — 175 MeV, and k = 0.18 is represented by the 
SETH column in Table 01. We see a good agreement with the data. The result of [is and k variations 



is also shown in Fig. 15 



1.25 



1.2 
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Figure 15: The predictions on the semileptonic branching ratio of B° meson and the average yield 
of charmed quarks in its decays with the ±3%-variation of hadronic scale in the nonleptonic effective 
lagrangian and the change of doubly charmed mode fraction k = 0.15 -j- 0.22 in comparison with the 
ARGUS and CLEO data shown by the dot with the error bars. 



We stress that the theoretical OPE prediction for the contribution of b — > ccs mode, k = 0.4, 
is significantly overestimated, to our opinion. Indeed, looking at the semileptonic decay with t + v t 
in the final state we see that at the same values of form factors as in modes with light leptons the 
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heavy lepton mode is suppressed because of the restricted phase space, so that the reduction factor 
equals 0.22. The same effect has to take place in the 6-quark decays with two charmed quarks in the 
final state. So, since the sum of D and K masses is greater than the r mass we could expect even 
the greater suppression than for the heavy lepton mode. Nevertheless, OPE operates with the quark 
masses at so moderate release of energy, that certainly leads to the overestimation of phase space in 
the decay of b — > ccs. We use more realistic values for k close to 0.2. 

Next, the contribution of four-quark operators are suppressed in decays of A b baryon. So, one 
expects that the deviation of its total width from that of B mesons is about 2-3% ||(|. To the 



moment, the experimental result is far away from this expectation (see Table ^). We can reproduce 
the data on the A b lifetime, if we introduce the light diquark scale 

fid = 4 A ' ( 5e ) 

where T/2 corresponds to the half tension of color string inside the diquark. Then we pose 

(57) 

which result in a good agreement with the experimental data. 

To finish this discussion we draw the conclusion: a probable way to reach the agreement between 
the theoretical predictions and available experimental data on the lifetimes and inclusive decay widths 
of heavy hadrons is to suggest the different normalization points in the effective nonleptonic lagrangian 
for the heavy quark weak decays as dependent on the hadron. This assumption provide us with quite 
acceptable results to the leading order in a s . The variation of normalization point shows the sensitivity 
of calculations to the higher orders in the QCD coupling constant, which indicates, first, the necessity 
to proceed with the higher orders, and, second, the appropriate choice of scale can allow us to decrease 
the scale-dependent higher orders terms. 

Thus, we suppose that the preferable choice of scale in the c — > s decays of B c is equal to 

/4=/^-/^-~(0.85GeV) 2 , (58) 
and at ai(/x.B c ) = 1.20 in the charmed quark decays we predict 

t[B c ] = 0.48 ± 0.05 ps. (59) 



7 Conclusion 

We have investigated the semileptonic decays of B c meson due to the weak decays of charmed quark 
in the framework of three-point sum rules in QCD. We have found the important role played by the 
Coulomb-like a s /v -corrections. As in the case of two-point sum rules, the form factors are about three 
times enhanced due to the Coulomb renormalization of quark-meson vertex for the heavy quarkonium 
B c . We have studied the dependence of form factors on the threshold energy, which determines 
the continuum region of bs system. The obtained dependence has the stability region, serving as 
the test of convergency for the sum rule method. The HQET two-point sum rules for the leptonic 
constant fs s and fs* have been reanalyzed to introduce the term caused by the product of quark and 
gluon condensates. This contribution essentially improves the stability of SR results for the leptonic 
constants of B mesons, yielding: /# = 140 -j- 170 MeV. 
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We have studied the soft limit for the form factors in combined HQET/NRQCD technique at the 
recoil momentum close to zero, which allows us to derive the generalized relations due to the spin 
symmetry of effective lagrangian. The relations are in a good agreement with the full QCD results, 
which means that the corrections to the form factors in both relative velocity of heavy quarks inside 
the be quarkonium and the inverse heavy quark masses are small within the accuracy of the method. 

Next, we have studied the nonleptonic decays, using the assumption on the factorization of the 
weak transition. The results on the widths and branching fractions for various decay modes of B c are 
collected in Tables. 

Finally, we have estimated the B c meson lifetime, and showed the dependence on the scale for the 
hadronic weak lagrangian in decays of charmed quark 

t[B c ] = 0.48 ± 0.05 ps. 

Our estimates are in a good agreement with the theoretical predictions for the lifetime in both the 
potential models and OPE as well as with the experimental data. 

This work was in part supported by the Russian Foundation for Basic Research, grants 99-02-16558 
and 00-15-96645. 

8 Appendix A 



For the perturbative spectral densities pi(si, s 2 , Q ) we have the following expressions |TB 



3 k 

p+(s l7 s 2 ,Q 2 ) = ^^{-(Ai + A 2 ) - /c[m 3 (m 3 - mi) +m 3 (m 3 - m 2 )] - 

[2(s 2 A 1 + s 1 A 2 )- M (A 1 + A 2 )] (60) 
■[ml - - + m 1 m 2 - m 2 m 3 - mim 3 ]}, 
3 k 

p-(si,s 2 ,Q 2 ) = -—^{-(Ai - A 2 ) - /c[m 3 (m 3 -mi) - m 3 (m 3 - m 2 )] + 

[2(s 2 Ai-siA 2 )+ M (Ai- A 2 )] (61) 

•[wij — — + m 1 m 2 — m 2 m 3 — mim 3 ]}, 
3 

Pv(s!,s 2 ,Q 2 ) = r^{(2siA 2 - uAi)(m 3 - m 2 ) 

+ (2s 2 Ai -MA 2 )(m 3 -mi) +m 3 A;}, (62) 

3 1 
Po(si,s 2 ,Q 2 ) = -^{( mi - m 3 )[ml + -( Sl A 2 2 + s 2 Aj - uA^)} 

-™2( m 3- y)-™i( m 3- y) (63) 
+m 3 [m3 - ^(Ai + A 2 - u) + mim 2 ]}, 

pi(s u s 2 ,Q 2 ) = ^^{m 1 [2s 2 A 1 -uA 2 + 4A 1 A 2 + 2A| 

mimg[4s 2 - 2u] + m 2 [2siA 2 - uAi] - m 3 [2(3s 2 Ai + s 1 A 2 ) 

-u(3A 2 + Ai) + k + 4AiA 2 + 2A^ + ml(As 2 - 2u)} (64) 

g 

+ -(mi - m 3 )[4sis 2 AiA 2 - w(2s 2 AiA 2 + siA^ + s 2 A\) 
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+2s 2 ( Sl Al + s 2 Al)}}, 
3 

P A (s 1 ,s 2 ,Q 2 ) = -^j^{kum 3 (2m 1 m 3 -2ml + u) + l2(mi-m 3 )slAj + 

kA2[(mi + m 3 )u — 2s\(m 2 — ^3)] + 2A\{k + 3usi)(m\ — m 3 ) 
^A^ku^ - m 3 ) + 2A 2 (k - 3u 2 ){m 1 - m 3 )] + (65) 
2s 2 (mi - m 3 )[2kml - jfeAj + 3uAj - 6wAiA 2 ] - 
2sis 2 {km 3 - SA^mi - m 3 ))], 

p^{s 1 ,s 2 ,Q 2 ) = -^{-2{m 1 -m 3 )[{k-3us 2 )A 2 1 + 6slA 2 ] + 

ku(mi — m 3 )(2m 3 + A 2 ) + ku 2 m 3 + A\\ku(2m\ — m 2 — 3m 3 ) 
-2(m 1 -m 3 )(ks 2 -kA 2 + 3u 2 A 2 )]- (66) 
2s 1 [(m 1 - m 3 )(2km 2 3 - QuA x A 2 - 3uA 2 2 ) + 
2s 2 {km 3 + 3m\A\ — 3m 3 A 2 ) + A;A 2 (2m! — m 2 — 3m 3 )]}, 

p' A (s u s 2 ,Q 2 ) = ^{2(m 1 -m 3 )[(k + 3us 2 )A 2 + 6s 2 Al) + 

ku{m\ — m 3 )(2m 3 + A 2 ) + ku 2 m 3 + Ai[ku(—2mi — m 2 + m 3 ) 
-2( mi -m 3 )(ks 2 -kA 2 + 3u 2 A 2 )}+ (67) 
2si[(mi - m 3 )(2£;m 3 - 6mAiA 2 + 3uA^) - 
2s 2 (km 3 — 3m 1 A 2 + 3m 3 A\) + kA 2 {2mi + m 2 — m 3 )]}. 

Here k = (si + s 2 + Q 2 ) 2 — 4s]S 2 , u = Si + s 2 + Q 2 , Ai = s\ — m\ + m 3 and A 2 = s 2 — ml + m\. 
mi,m 2 and m 3 are the masses of quark flavours relevant to the various decays, see prescriptions in 
Fig. |. 

9 Appendix B 

Here we list the expression for the form factors of semileptonic decays B c — > B^*> taken from the 



potential model [32 . 



f = _ ^-^;^) ^w, m 

F ° = 2Al Vm^ W ' (70) 



= U2,i. t /M t jfe. 

2m s V Mb '• K ' 
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where 

2loco x / 2ljl) x I ml(w 



eW = Z^fV^T^ exp p^T^J' (73) 

^ = 2^^-j , ^. = 2^ ^ j , (74) 

where w is the product of B c and .B^ four-velocities. The quark masses and the leptonic constants 
have the values generally used in the calculations in the framework of potential models 

m b = 4.8 GeV, m c = 1.5 GeV, m s = 0.55 GeV, f Bc = 0.47 GeV, L w = 0.17 GeV. 
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